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Abstract 

We introduce a class of stochastic weakly coupled map lattices, as models for studying heat conduction 
in solids. Each particle on the lattice evolves according to an internal dynamics that depends on its 
energy, and exchanges energy with its neighbors at a rate that depends on its internal state. We study 
energy fluctuations at equilibrium in a diffusive scaling. In some cases, we derive the hydrodynamic 
limit of the fluctuation field. 
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1 Introduction 



It is generally expected that the temperature profile of thermically isolated solids relaxes according to 
some non-linear heat equation 

dtT = v(L»(r)vr). 

The derivation of this equation, starting from a microscopic hamiltonian dynamics, is one of the major 
open problems in statistical mechanics out of equilibrium [5] |14) . Very idealized models of solid constitute 
a possible starting point to develop our understanding on the problem, and in fact, at the present time, 
even such cases can be challenging to analyze. 

Lattice gases furnish the first examples where the heat equation can be recovered through a diffusive 
rescaling of space and time. Each atom reduces here to a point on a lattice, is characterized by its 
energy only, and interacts stochastically with its neighbors. The stochastic interaction is a reasonable 
approximation for particles that should evolve according to some chaotic internal dynamics. A set of 
independent random walks (IRW), the simple symmetric exclusion process (SSEP) [10,, and the Kipnis 
Marchioro Presutti model (KMP) 11 , constitute three cases where the heat equation can be easily 
derived. These systems may be closer to physics than what could be thought at the first glance: in [12) . 
an example of noisy hamiltonian system is found, that gives rise to the SSEP in a weak coupling limit. 

This said, true atoms should clearly be described by position and momentum, and not only energy. 
To move towards a more realistic situation, let us thus consider a lattice where particles at the nodes 
now have some internal degree of freedom. We assume the local dynamics of particles to have good 
ergodic properties, and the interaction between them to be weak and controlled by a parameter e > 0. 
So, individual atoms will reach their own equilibrium at a smaller time scale than energy is exchanged. 

Much progress has been recently accomplished in the understanding of these systems, by means of a 
two-stage analysis, as outlined in [S]. First, letting e and rescaling time properly, an autonomous 
stochastic process for the energy alone is obtained. Second, rescaling space and time diffusively, the heat 
equation is derived starting from this mesoscopic stochastic process. The first part of this program has 
now been rigorously accomplished for a purely hamiltonian dynamics [5 . Though in a rather different 
way, the distinction of two different time-scales also has allowed to better understand the conductivity 
of weakly anharmonic solids [T] [T3] . 

The introduction of two successive limits is however somehow artificial. One eventually would like 
to fix some, possibly small, e > 0, and directly condsider the diffusive limit. Results in that direction 
are few. Let us quote three of them. First, in [2], the heat equation is recovered starting from a noisy 
harmonic chain (see also |7l for fluctuations at equilibrium). The proof rests on the fact that an exact 
fluctuation-dissipation decomposition is available there. Second, a deterministic coupled map lattice is 
analyzed in [3] , and diffusion of an initially localized energy packet is shown. The dynamics of particles is 
there assumed to be independent of energy. Third, fluctuations of energy for a noisy hamiltonian system 
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at equilibrium are considered in [T5] (see also [5] for a simpler model) ; the flucuation field is shown to 
converge to a generalyzed Ornstein-Uhlenbeck process. Some level of symmetry of the generator, the 
sector condition, is needed to derive this result. 

In this paper, we introduce models of heat conduction that are conceptually very simple. They maybe 
will even look a bit cheap compared to the aforementioned examples. Though, the analysis reveals that 
nowdays techniques do not directly apply to describe how energy diffuses in the cases we consider. We 
will in fact need to look at particular situations in order to derive rigorous results. To make this more 
concrete, let us now briefly outline the content. 

In Section [51 we define three models of conduction. In a sense made precise by (|2.3p below, they 
constitute respective versions of the IRW, SSEP and KMP models, where, now, atoms have an internal 
degree of freedom. The interaction between near particles is controlled by a parameter e > 0. The models 
so obtained are non-gradient and do not satisfy the sector condition. 

Like in [IS] and [S] , we study the fiuctuations of energy in equilibirum, as this might be the easiest 
problem to look at. We introduce the fluctuation field in Section |31 and we state there Theorem [1] 
that will allow us to derive the hydrodynamic limit of this field in some cases. While it is essentially a 
variation on a classical method 10 , we here assume that exists an approximative fluctuation-dissipation 
decompositon of the microscopic current in a particular sense. We hope that this kind of result can be 
of some use elsewhere. The proof is postponed to Section \5\ 

Explicit results are gathered in Section 2] The hydrodynamic limit of the field in the first order in e is 
derived in Theorem [21 As a corollary, a weak coupling limit is shown, where e is sent to after a diffusive 
rescaling of space and time. This kind of limit is a bit more satisfactory than the one in [5] [12] [S] , where 
e was sent to for a fixed number of particles without any rescaling of space. It also avoids the two-stage 
procedure of [Sj. One then approximates, like in 4 , the local dynamics by a process independent of 
energy. In Theorem [31 the hydrodynamic limit of the flucuation field is obtained under that hypothesis 
for any small enough e > 0. While, with respect to [4 , we deal with a rather particular case, our proof 
is much shorter, and we are not limited to an initially localized energy packet. 

2 Models 

Let us consider a set of > 3 identical particles. The coordinates of the particles are indexed by a point 
in the periodic one-dimensional lattice Z^r of integers modulo N. They form thus a one-dimensional 
chain with periodic boundary conditions. Each particle is characterized by a local degree of freedom, and 
a non negative energy, so that the phase space X is given by 

X = X M^. 
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where 



M = Tx{-l,l} = (R/Z) X {-1,1}. 
A point r G M is written as r = (g, i;) e T x { — 1, 1} ; a point x e X is written as 

X = (r, e) = (ri, . . . ,rjv,ei, . . . ,ejv) = (q,^, e) = (qi, . . . , gAr, ^i, • . • , ^at, ei, . . . , ejv)- 

We set also Xk = {rk,ek) = {qk,^k,ek) for 1 < fc < A^. 

We see qk as the position of particle fc, and <;k as the sign of its velocity. It could have seemed more 
natural to first define the usual hamiltonian variables qk and pk , and then express the energy in function 
of them. However, for our purposes, the distinction between the local variable and the energy Ck is 
the most important. In fact, we almost never will make use of the variables qk and ^fc. We assume that 
each particle have a ground state energy eg > 0, and we think of Ck as being the extra amount of the 
energy of particle k with respect to Bq. 

Let e > 0. We consider a stochastic dynamics defined through the generator 



acting on functions on the phase space X. Here L^c governs the internal dynamics of particle fc, and L'l. 
the exchanges of energy between particles k and fc + 1. We denote by Xt the process generated by L at 
time t > 0. Let us now concretly define L' and L" . 

2.1 Local dynamics 

Let 1 < fc < iV. One defines 

L'f.u = A'f.u + S'kU = qk Veo + Ck-dq^u + Vso + Sk ■ {u{. . . , -<;fe, . . . ) - u(. . . , <jfc, . . . )) . 

Here AJ, and S*^ denote respectively the symmetric and antisymmetric part of L'^,, with respect to the 
uniform probability measure /iu on M. So particle fc moves freely on the circle at speed \/eo + Ck and, 
at random times, the sign of its velocity gets fliped. We will make use of the following notation: given a 
function u on X, one writes 



No sector condition holds. Moreover, the dynamics generated by the symmetric part 5^ alone is very 
degenerated since it does not affect the physical location qk- Up to irrelevant details, the convergence to 
equilibrium of the dynamics generated by L'j, on M has been studied in [12] , where an explicit expression 
for the transition probability is given (see also [6] for an approach based on hypocoercivity). One has 
the following for 1 < p < +oo. There exist constants c > and C < +oo such that, for any u E LP^(X) 
satisfying (u|e) — for all e, and for any t > 0, 



L = L' + eL" = {L'k + eL'D 




||e^'*u||p < Ce-=*||^.|| 
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We will use this information in the following form. Let / € Lp^^(X) be such that (/|e) = for all e. 
Then there exists a unique function u E L^^(X) satisfying {u\e) = for all e, that solves the Poisson 
equation 

-L'u = f. (2.1) 

We call this function u the fundamental solution to this Poisson equation. Moreover there exists a 
constant C < +oo independent of / such that 

\W\\p < C||/||p. (2.2) 

This follows from the explicit representation 

Jo 

2.2 Interactions 

In a one-dimensional chain of oscillators evolving accoring to the laws of classical mechanics, the strength 
of the force between two atoms depends on their relative physical location, and the instantaneous current 
between them is a function their positions in phase space. To mimick this feature in our model, let us 
introduce a uniformly bounded symmetric function x > on M^. The function x 'v^^iU niake the rate of 
exchange of energy depend on the locations of particles. We set also 



X = / x(^o,?-i)Mu(dr). 

As an example, one can take 

x{rQ,ri) = x{qo,qi) = XlaM^ilo,qi), 

where X[a,b] denotes the characteristic function of some interval [a, b] C T. Then, the interaction between 
two particles will only be possible if they both are in the interval [a, b]. 

We will define three possible ways of interacting, and so three different models. In the three cases, 
the interaction will be purely stochastic, meaning that the operator L" is symmetric with respect to the 
canonical invariant measures defined in Section [31 This is admittedly in sharp contrast with hamiltonian 
interactions. We may view our models as perturbations of the IRW, SSEP or KMP dynamics respectively. 
Indeed, the operator C acting on functions u that depend only on the energies e, and defined by 

Cu = ^ {L"u\e), (2.3) 
X 

will be nothing else than the generator of these dynamics. 
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Model 1: Exchanging bits of energy. Let us restrict the set of reachable energies toN = {0,l,2,...}, 
so that the phase space writes X = x N^. Let 1 < k < N. One defines 

L'kU = x{rk,rk+i) ■ (efe(u o - u) + ek+i{u o cr^^^^ - u)) 

with 

<^fe W = (r,ei, . . . ,efc_i + l,efc - 1, e/s+i, . . . , ejv), 
^ki^) = . • . ,efc_i,efe - l,efe+i + 1, . . . ,ejv)- 

So, with such an interaction, particles tend to give part of their energy to their neighbors at a rate 
proportional to the value of their energy times the value of x- The instantaneous current is defined 
through the relation 

Lck = eL"ek = e(jk-i,k - jk,k+i) 

with 

jk,k+i = x{rk,rk+i){ek - Gk+i)- (2.4) 

Model 2: Exchanging the whole energy. Let us further restrict the set of reachable energies to 
{0, 1}, so that the phase space now writes X = x {0, 1}^. Let 1 < k < N. One defines 

L'kU = x{rk,rk+i){uo ak - u) 

with 

UfeX = (r, ei, . . . ,efe+i, e/s, . . . ,ejv)- 
So here, near particles may exchange their energy. It is checked that the instantaneous current is the 
same as in the first model. 

Let us stress at this point that the dynamics still should perfectly make sense if energy was allowed 
to take more than two values. However, since energy is only exchanged between particles, there should 
be n — 1 conserved quantities if energy was allowed to take n values. It is not at all clear in that case 
that an initial density of energy should evolve autonomously in the hydrodynamic regime, since an initial 
macroscopic state should in fact be specified by the value of each conserved quantity. We have thus 
limited ourselves to the case of two reachable energies only, in order to avoid extra complications. 

Model 3: Sharing energy. Energy is here allowed to take any non negative value, so that the phase 
space is given by X = x M^. Let 1 < A; < iV. One defines 

L'fe'u = 2x{rk,rk+i){Tk-Id)u 

with 

TkU = / u(r, ei, . . . ,p(efc + efe+i), (1 -p)(efe + efc+i), . . . ,eAr)dp. 
The instantaneous current is the same as in the two previous models. 
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3 Fluctuation field 



The derivation of the heat equation requires to consider an energy profile initially out of equilibrium. In 
this work, like in [TS] [S] , we will however only look at small perturbations out of equilibrium. To see what 
small means, let us assume that our system is in equilibrium at an inverse temperature < /3 < +00 
(see below), and let us denote by (e)^ the average energy of a particle. Let / C T be an interval. For a 
typical event in the invariant measure, one has 

(e..-(e),) ^ Vn. 

k:k/N£l 

So let 

k:k/NeI 

be the macroscopic observable we look at. The random variable £t{I) is expected to fluctuate in a diffusive 
time-scale. So a typical kind of probability that one would like to estimate in the limit — > 00 is 

P(£jv=,(/)~6|£o (/)-«), 

where a, 6 G R, and where a; ~ a means x G [a — S, a + S] for some small 6 > 0. So in fact, we are 
considering intial profiles that have an excess of energy of order ^/N in subintervals of T, rather than of 
order N for a profile truly out of equilibrium 

We now need some definitions. Let us first define canonical equilibrium measures. For our three 
models 

{u)f, ^ / u(r,e)/iu(dr)/x^(de), 
Jx 

where, for u = u{e), one has 

I_il^{u) = e-'3('=i+-+*=") u(e) for Model 2, 

^p{u) = ^— / e-^('=i+-+^")w(e)de for Model 3, 

with Z[j3) a normalization factor. 

Let us next define the fiuctuation field (Y^^^) ^^^^^j,. For any smooth test function H on T, we define 

1 ^ 

rW(^) ^ --YH{k/N){e,{Nh)-{e)p). 

Most of the time we will just write Y instead of F'^'. Formally, {YN^t)o<t<T is as a random process 
with values in I?([0, T], "H-a), where T) denotes the set of cad-lag functions, and where H-3 denotes the 
Sobolev space of distributions Z on T with norm ||Z||^3 = X^fcez with ek{x) = e^"*^^. 
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Let us finally define the stationary generalized Ornstein-Uhlenbeck process (Yt)o<t<T € C([0, T] , 
where C denotes the set of continuous functions, for some diffusion constant eD > and thermal capacity 
(T^ > 0. We here have anticipated the fact that, in our applications, the diffusion constant will be of 
order e. Let H he a, smooth test function on T, and, for t > 0, let Ht be the solution to the heat equation 

Ho^H, dtHt = tD.dlHt. (3.1) 

First, Yo(-ff) is a centered Gaussian random variable with variance 

cr^ / Hl{x)Ax. (3.2) 

Next, for t > 0, 

Yt{H) = Yo{Ht)+Mt{H) 
where Mt{H) is a Gaussian centered martingale with variance 

2€Da^ {d.^H,{x)f dx^ ds. 

Our aim is to show that {Yj^2f)t converges in some sense to {Yt)t as N ^ 00, with 

D = D{P) > and = a^{P) = {3^)13- (e)}. 

We say that a function D . Ve on X is gradient if it is of the form 

D.Ve^ ^ Z?fc.(efc-efc+i). (3.3) 

We then also write D = X^fcGZjv ^'^^ This quantity depends on the representation p.3p of the gradient 
function D . Ve, which is not unique due to periodic boundary conditions. This however will create no 
trouble in the sequel. We will use the following theorem, shown in Section [Sj that furnishes a possibly 
interesting bound when e' << e. The hypotheses are close to optimal for the proof, but could be quite 
simplified for our applications. 

Theorem 1. Let e' > 0. Let /3 > 0. Let a'^{P) = (e^) 

/3 ~ (^)|- Assume that there are sequences of 
functions ((_D . Ve)Ar)j^, (ujv)a' and {gN)N on X, that may depend on j3, such that 

-LuN - ^{jo.i ~ {D .Ve)N) = e'L"gN, 

that 

J™ Tn \{TkUN,UN}f3\ = 0, 

N—^oo iV — ' 

fceZjv 

and that the following quantities converge as N 00: 
E K^^)^l' E \{9N,r,L"g^),l E 

\{uN,jk,k+l)l3\, E \{''^N,TkL"gN)p\- 
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Let then 



N^oo ^ — ' 

Let also H be a smooth test function on T. The fixed-time distributions of the process {Ypf2f(H))t converge 
to the fixed-time distributions of a process (Yt{H))t given by 

Yt{H) = Yo{Ht) + Mt{H) + Rt{H), 

where Ht is given by p.ip . where Yo(i?t) is a Gaussian random variable with variance given by p.2p 
with Ht in place of Hq, where Mt{H) is a Gaussian martingale which variance is given by 

2e / ( f{d,H,fdx]ds.(xa^il3)~ lira V (u^, Jfe,fe+i)^i + - lim {uN,TkL"gN)p) , (3.4) 
and where Rt{H) is centered with variance bounded by a constant times 

+ lim J2 \{9N,nL"gM)p\. 

4 Diffusion of energy 

The results below hold for our three models. We first obtain the hydrodynamic limit of the field Yjqi^ to 
the first order in e, with a bound on the variance of the rest term. 

Theorem 2. Let D — x- Let also H be a smooth test function on T. The fixed-time distributions of the 
process (Ypf2f.(H))t converge to the fixed-time distributions of a process {Yt{H))t given by 

YtiH) = YoiHt) + Mt{H) + RtiH), 

where Ht is given by p. II) . where YQ{Ht) is a Gaussian random variable with variance given by p.2[) 
with Ht in place of Hq, where Mt{H) is a Gaussian martingale which variance is given by 

2eDa^P)J^ (^J{dM^dx^ds + 0{eH), 

and where Rt{H) is centered with a variance that is 0(6^(1 + t^)). 

Proof. Thanks to (|2.ip . there exists a bounded function v depending only on xq and xi, solving the 
Poisson equation 

-L'v ^ jos ~ D . {eo ~ ei) = (x - x)(^o, ^'i) ■ (eo - ei), 

and such that {v\e) = for every e. Defining e' — e^, un — and gjv = —v, one sees that the hypotheses 
of Theorem [T] are satisfied, from which our claim follows. □ 

Thanks to the good bound in 0{£^{1 + on the variance of the rest term, a weak coupling limit 
is obtained as a corollary of Theorem [51 Let {Yt)t be the limit process appearing there, where we now 
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assume that < i < e~^T. Since T is arbitrary, there is no difficulty in doing so. Then, as e — 0, the 
fixed-time distributions of the field (Y^-if.)o<t<T converge to the fixed-time distributions of a generalized 
Ornstein-Uhlenbeck process with diffusion constant D = x and thermal capacity a'^(j3). 

Like in [4j , we next approximate the internal dynamics of particles by a process independent of energy. 
So we now assume that, for fc G Zjv, 

L'kU = i,kV^-dq^u + Vsq. (-u(. . . , -<^fe, . . . ) -?/(..., <rfe,...)). (4.1) 

This is expected to be a reasonable approximation as long as energy takes values not too far from Sq. 
With such a definition of L', particles evolve now completely independently from each others. This in 
fact is not needed ; any operator i', not involving energy and generating a dynamics that quickly relaxes 
to equilibrium independently of the system size N , should be fine for our purposes. 

Theorem 3. Assume that L'/^ is given by (14. ip for k G Z^v, and that e > is small enough. The fixed- 
time distributions of the field {Ypf2f.)t converge to the fixed-time distributions of a generalized Ornstein- 
Ulhenbeck process with some diffusion constant eD > independent of the temperature, and with the 
thermal capacity = . 

Proof. The proof is made of two steps: we first show that Theorem [T] can be applied with e' = 0, and we 
next identify the variance of Mt{H). 

Step 1. Let n>l. We show that the hypotheses of Theorem [T] are satisfied with e' = e"+^, and that 
the limit as — > cx) of all the quantities appearing in the hypotheses of Theorem [TJ are bounded by a 
constant that does not depend on n. Let us introduce an extra definition. We say that a function u on 
X is pre-gradient if it is of the form 

M(r,e) ^ ^ u(r, fc)(efc - Cfc+i). 

Let us first make two observations. Given a pre-gradient function /, the fundamental solution u to 
the equation —L'u = / is itself pregradient. Indeed the coefficients u{y, k) are the fundamental solutions 
to the equation —L'u{r,k) — /(r, fc). Second, if a function u is pregradient, so is L"u. Indeed, for the 
three models, one has 

L"u{r,e) = {u{r,k-l)-2u{r,k)+u{r,k + l))x{rk,rk+i){ek-ek+i). (4.2) 

So moreover, if u depends on variables x^, . . . , a^b with a <h, then L"u depends on variables Xa~i, . ■ ■ , Xh^i. 

Now, for N large enough for given n, the functions un, {D . Ve)^ and g^r that we will construct will 
actually not depend on A^, so that we drop the index TV. We set 

n n 

u = Ye''u^^\ (D.Ve) = ^e^-^ D^''\Ve and g = (4.3) 

k=l k=l 
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with the following recursive definitions of u'*^) and D^'^) for every fc > 1. First 

I?(i).Ve = Oo,i|e) =x(eo-ei), 
and u^^^ is the fundamental solution of 

-L'm^^) ^ jns - Oo,i|e)- 
The function u*^^-* is pregradient. Next, given the pregradient function ^i*^*^^ for fc > 1, one takes 

^(fc+i).Ve= (^''TiWle)^ 

This indeed defines a linear gradient since L"u'^^^ is pregradient. Then one defines u'^^^^^ to be the 
fundamental solution to 

The function is pregradient. 

The functions m, D . Ve and g defined by (|4.3p satisfy 

-Lu~ e{jo,i ~ D . Ve) - e"+^L"g. (4.4) 

Thanks to (|2.ip and the convergence of geometric series, the quantities appearing in the hypotheses of 
Theorem [1] are then indeed bounded independently of n if e > is taken small enough. 

Step 2. Let u, D . Ve and g be given by (j4.3p Let us first give an expression for D. Since {u\e) — 
for every e, one has also {L'u\e) — for every e, so that, by (|4.4p . 

D.We ^ x{eo-ei) + {L"u\e)+e"{L"g\e). 

The function L"u is computed by means of (|4.2p . Writing 13 . Ve = J2jei.N ^i^^i ~ ^j+i)' '^'^'^ 

D = ^1 = X+ E (x(-)(«(-,fc-l)-2u(-,fc)+«(.,fc + l)))+0(e"). 

Taking on the other hand the expression p.4p for the variance of Mt{H), and using (|4.2p to compute 
limAT^oo Z^fceZjvi^JV, jfc,fc+i), one gets 

a'^D = a^x- lim ^ (u^, + 0(e"). 

N~>oo ^ — ' 

This shows our claim. □ 

5 Proof of Theorem [1] 

We first remind a general formula, that we will use several times. If u = u{x,t) is a regular enough 
function on X x M+, then the process 

Mt = f Lu(X„s)ds+ f dtu{Xs,s)ds + uiXo,0)-uiXt,t) 
Jo Jo 
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is a centered martingale, which variance at equihbrium is given by 

pt 

E^{MtY = 2 / {ui;s),i~L)u{;s))pds. 







Let us then make some comments on the notations. We will write (•) for For k G Z^r, we define the 
centered variables 

rjk = efc - (e). 

We will write gradient functions as D . V77 instead of D . Ve. We will not write explicitly the dependence 
in N of the functions un, gN and (D . appearing in the hypotheses of Theorem [1] We will use the 

same symbol dxf for both the true derivative of / and its apporximative derivative N (^f {x + 1 / N) — f (x)). 
Proof of Theorem [H Let H = H{x) be the smooth test function given in the hypotheses, and let 
G = Gix, s) be defined for < s < NH, by 

iV^StG + . = 0, G{-,NH) = H{-). (5.1) 
Since Lr]k = e{jk~i,k ~ jk,k+i) for k G Zjv, one has 

/ E G{k/N,s)L^kds = ^ E G{k/N,s)e{jk-i,k~Jk,k+i)ds 

There exists thus a martingale M^l^ such that, after an integration by parts, 

YN2,{G{;Nh))-YoiG{;0)) = -a4%+ / — V dtG{k/N,s)r,kds 



1 /"^ * 1 

iv / ~7W E <9a:G(fc/iV, s)ejfc,fc+i ds. 



k£Zi, 

Using then the approximate fluctuation-dissipation decomposition of the current stated in the hypotheses 
of Theorem [1] one gets 

Y^2,{Gi-,NH))-Yo{Gi-,0)) - -a4%+ / -j= V dtG{k/N,s)rjkds (5.2) 

■^0 VN 



Jj [ ^ E d.Gik/N,s)Tkuds 

N 7^ E ^-G(fc/iV, s)r,I? . Vr; ds 



(5.3) 
(5.4) 
(5.5) 



We now proceed to the analysis of the last three terms in the right hand side of this equation. Let 
us start with (j5.3[) . Since dtdxH{k/N, s) = — (eD / N'^)dxH {k / N , s), there exists a martingale M^l^ such 
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that 



The hypothesis that X^fceZw '")j9l remains bounded as — > oo, ensures that E(pj^2j)^ — > as 

TV — )■ oo. Let us then look at (j5.4p . Reminding that Z? = hmAr_i.oo X^fcGZjv integration by parts 
yields 

eD f^'* 

(EH) =]^y^ r,(5^G(.,5))ds 



eD 
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='■^1 n(5^ij(.,s))ds + p«,. 

Using now that X^feeZw l^*^! remains bounded as iV — oo, and that 

92G(fc/7V,s)-a2G((fc+j + l)/iV,,s) = d',G{k/N,s)+0{\j + l\/N), 

one deduces that E(pj^2j)^ — as — oo. 

We finally handle the less usual term (j5.5|) . which is a centered process. It is convenient to define 

= E d.G{k/N,s)T,g, ^ d^Gik/N, s) 

For z > one then defines Vz,s and V^^s as repsective unique solutions to the equations 

(z-L)K,. - -L"F,, {z-L)V,.s = (5.6) 
So, defining a centered martingale M^a^, one has 



1 ' 

(ESI) / (-L")Gs(X,)ds 

Therefore, taking z — 1/iV^, one finds a constant C < +oo such that 
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But one has the bounds 

z{Vl:) < e-\FsA-L")Fs), {V.,sA-m,s) < e-\F,,{-L")F,) , (5.7) 
as well as the same bounds with the tilted quantities. Indeed, (j5.6l) implies 

< 6-1/2 i^p^ i^L")F,y/^V..s . {-eL")V,.,Y'^ 

< £-1/2 (F, . (-L")F,)V2(v;_, . (-L)K^,)i/2, 

from where (|5.7p is derived. So finally we obtain 

^pim') < C{ere-\l + f )\im ^ \{g,L"ng)p\. (5.8) 

So at this point, let us define 



N- 



RnH = (Ell) and MN2t = -M^X-^'^NH- 
Taking into account that G satisfies (|5.ip . one has so far obtained 

Since E^(i?^2j) is bounded by (|5.8p and since both p^a^ and p^2( converge to in L'^{Pp), the proof will 
be concluded by resacling time, if one establishes that, as — > oo, 



Ep{M^,t) ^ 2e^lim^^ f f f {d^Gfix,s)dx\ds.(x<^^ - ^ {u,jk^k+i) + ^ 

Jo \J / fcezjv ^ fceZw 

Wrting for Kj(^f(-,s)) and putting 

= jf^T. d.Hik/N,s)TkU, 

the variance of Mjqif is given by 

E(Mjv2t)' = 2 / + C/„ (-L)(r. + Us)) ds 

Jo 

l-N^t i-N^t i-N^t 

= 2e {Ys,{-L")Ys)ds + 4e {Us, {-L")Ys) ds + 2 (C/,, (-L)C/,) ds. 
Jo Jo Jo 

One first readily computes that, since X^feeZjv li'^^ jk,k+i) \ remains bounded as — >■ oo, 
2eJ^ {Ys, i-L")Ys) ds ^2ex<y^ . lirn^^ (^Jid,G)^{x,s)dx^ds, 
46 / {Us,{~L")Ys)ds ^ -4ejini^— / ( id,G)^{x, s) dx\ ds . i^^w)- 

•'^ Jo \J J 



(5.9) 
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Then, one uses the approxiamte fluctuation-dissipation decomposition of the current stated as hypothesis 
in Theorem [1] to obtain 

k^I^N k^ljf-i kG^N 

The second term in the right hand side of this equation wiU not contribute since it is a gradient. Tak- 
ing into account that both X]fegz„ ^^'^ SfegZjv remain bounded as iV -> oo, A 
computation furnishes 

2 / {Us, i-L)Us)ds ^ 2 Jirn^— / ( / (d^Gfix, s)dx\ ds . (e {u, Jk.k+i)+e' J2 i'^^^kL"g) 
So ([F^ is shown. □ 
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